Band mass anisotropy and the intrinsic metric of fractional quantum Hall systems 
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It was recently pointed out that topological liquid phases arising in the fractional quantum Hall effect (FQHE) 
are not required to be rotationally invariant, as most variational wavefunctions proposed to date have been. 
Instead, they possess a geometric degree of freedom corresponding to a shear deformation that acts like an 
intrinsic metric. We apply this idea to a system with an anisotropic band mass, as is intrinsically the case 
in many-valley semiconductors such as AlAs and Si, or in isotropic systems like GaAs in the presence of a 
tilted magnetic field, which breaks the rotational invariance. We perform exact diagonalization calculations 
with periodic boundary conditions (torus geometry) for various filling fractions in the lowest, first and second 
Landau levels. In the lowest Landau level, we demonstrate that FQHE states generally survive the breakdown 
of rotational invariance by moderate values of the band mass anisotropy. At 1/3 filling, we generate a variational 
family of Laughlin wavefunctions parametrized by the metric degree of freedom. We show that the intrinsic 
metric of the Laughlin state adjusts as the band mass anisotropy or the dielectric tensor are varied, while the 
phase remains robust. In the n = 1 Landau level, mass anisotropy drives transitions between incompressible 
liquids and compressible states with charge density wave ordering. In n > 2 Landau levels, mass anisotropy 
selects and enhances stripe ordering with compatible wave vectors at partial 1/3 and 1/2 fillings. 

PACS numbers: 73.43.Cd, 73.21.Fg, 71.10.Pm 



I. INTRODUCTION 

Two-dimensional electron systems (2DES) placed in a high 
magnetic field exhibit a wide variety of strongly correlated 
phases, which have been the subject of numerous theoretical 
and experimental investigations since the first observation of 
fractionally quantized Hall conductivity^. Examples of such 
phases of matter are the Laughlin states 2 ., describing partial 
fillings ^ = 1/3, 1/5 of the lowest (n = 0) Landau level (LL), 
as well as their generalizations to other odd-denominator fill- 
ings in the framework of hierarchy- and composite fermion 
theory 4 -. These phases are topologically ordered and possess 
quasiparticle excitations with fractional statistics. At half fill- 
ing of the first excited, n = 1 LL, an even more exotic paired 
state, the Moore-Read Pfaffian 5 , might be realized, which pos- 
sesses non-Abelian excitations - the Majorana fermions^. 

Besides the incompressible liquids, some fillings v also 
lead to compressible phases without quantized conductance. 
This is the case with the simplest of all fractions - v = 1/2 in 
n = LL - which is a Fermi liquid of composite fermions 7 , 
that only supports an anomalous Hall effect. Generically for 
any v, apart from the incompressible liquids, the natural can- 
didates are compressible phases that break translational sym- 
metry, such as charge density waves (CDWs) 9 . Those were 
in fact proposed to describe the ground state of 2DES before 
FQHE was observed^. When v is very small (under 1/7), a 
correlated Wigner crystal indeed becomes energetically supe- 
rior to a Laughlin-type stat o 11 - 12 . Furthermore, when v > 2, 
several varieties of states with broken translational symmetry 
become energetically favorable. Around half filling of n > 2 
LL, the ground state becomes a charge density wave in one 
spatial direction or a "stripe"—^; away from half filling, two- 
dimesional crystalline order sets in, resulting in a "bubble" 
phase d 15 ' 17 . Some of these phases also occur in n < 2 LLs 



when the hard-core component of the effective interaction is 
significantly softer than Coulomb 18 . More recently, the ex- 
periments have shown 19 that it is possible to have, at the same 
time, the quantization of resistance and anisotropic transport, 
suggesting a possible coexistence of an incompressible liquid 
with a compressible ("nematic") phased. 

Theoretical understanding of the FQHE was pioneered by 
Laughlin's method of many-body trial wavefunctions^. Model 
wavefunctions can be formulated using the conformal field 
theory^, and conveniently evaluated in finite-size systems via 
exact diagonalization of the parent HamiltoniansSi. In ad- 
dition, excitation spectra containing quasiparticles/quasiholes 
above the ground state can be studied. Such analytical and nu- 
merical studies are made much easier by exploiting symmetry 
and the corresponding quantum numbers to characterize the 
ground state and excitations. To this end, rotational symmetry 
has been very useful 2 ^; ocassionally, periodic boundary con- 
ditions have also been use d 23,24 . 

However, as it was recently pointed out 32 , rotational sym- 
metry is not fundamental to the appearance of FQHE. In a the- 
oretical treatment of the FQHE, it is important to distinguish 
several "metrics" that naturally arise in the problem. The band 
mass tensor yields a metric that defines the shape of the LL 
orbitals. A second metric derives from the dielectric tensor 
of the semiconducting material, and defines the shape of the 
equipotential contours around an electron. Rotational invari- 
ance means that these two metrics are congruent, however in 
a real sample they might be different from one another, thus 
lifting the rotational invariance. 

It turns out, however, that a given FQH state also possesses 
an intrinsic metric that is derived from the two types intro- 
duced above. FQH fluids can be described as condensates 
of composite bosons^, which are topological objects that ex- 
plain the quantization of the Hall conductance and the emer- 
gence of fractionally-charged quasiparticles. However, apart 
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from topology, composite bosons also have a geometrical de- 
gree of freedom - the intrinsic metric that controls their shape. 
In a rotationally invariant case, the intrinsic metric is equal to 
the metric in the Hamiltonian; more generally, as we explicitly 
demonstrate below, the shape of composite bosons can be de- 
fined even in systems without rotational invariance. The fluc- 
tuation of the intrinsic metric plays an important role for the 
geometrical field theory of the FQHEs£, and determines the 
energetics of quasiparticles, collective modes etc. Generaliza- 
tions of the Laughlin wavefunction to the broken-rotational- 
symmetry case have been proposed for liquid crystal and ne- 
matic Hall phases- 7 , and very recently Laughlin and Moore- 
Read wavefunctions (as well as their parent Hamiltonians) 
have been formulated for the anisotropic case^. 

The motivation for studying the effect of anisotropy in 
FQHE is twofold. On the one hand, anisotropy probes the 
variations of the intrinsic metric of FQH fluids, a fundamen- 
tal physical quantity that relates to the geometric description 
of FQHE. Secondly, we explore the possible effects resulting 
from tuning the rotational-symmetry breaking by an external 
parameter. Note that the rotational invariance is explicitly bro- 
ken in real samples due to the presence of impurities, which 
are essential for the emergence of FQH plateaus. Further- 
more, it is possible to induce the breaking of rotational invari- 
ance by tilting the magnetic field^, or by using systems with 
anisotropic bands, e.g. many-valley semiconductors like AlAs 
or Si in the presence of uniaxial stress. The former method 
is performed routinely and belongs to the most popular tech- 
niques for studying the FQHE; the latter method is relevant to 
AlAs^ and some new classes of materials where FQHE might 
be studied. We provide brief arguments how the tilting of the 
field can be mapped to an effective variation of the metric, and 
then focus on the second method. 

This paper is organized as follows. In Sec. |II]we introduce 
and motivate the model for a FQH system with band mass 
anisotropy. In Sec. [HI] we discuss the intrinsic metric of the 
Laughlin state. We define a family of the Laughlin wavefunc- 
tions characterized by the varying shape of their elementary 
droplets. Intrinsic metric is determined variationally by opti- 
mizing the overlap between this family of wavefunctions and 
the exact Coulomb ground state. In SeclIVIwe perform exact 
diagonalization of finite systems at several filling factors to 
explore quantum phase transitions that occur as a function of 
the anisotropy. Our conclusions are presented in SecfVl 

H. MODEL 

Consider an electron moving in the plane with the perpen- 
dicular magnetic field Bz = V x A(r). The Hamiltonian can 
be written in the following, manifestly covariant form: 

K a <* = ^-g ab 7Ta7T b . (1) 

2m 

Here 7r a = p a — -A a (r),a = x,y represents the dynam- 
ical momentum, and g is the mass tensor parametrized by 
the anisotropy a and the angle of the principal axis cf>. The 
mass tensor is unimodular dct g = 1. In the isotropic case 



when g is the unit matrix, we can obtain the single-particle 
energies (Landau levels) by choosing, for example, a sym- 
metric gauge A x = By/2,A y = ~Bx/2. In this case, 
the dynamical momenta become tt x = —ifo-^ + ^-y and 

Tty = ~ifo-§^j — m terms °f tne magnetic length £b = 

\JhjeB. The Hamiltonian can be transformed into diagonal 
form K = (a 1 a + |) with the help of ladder operators 
a oc ir x + iiiy and a) oc ir x — iir y . However, for each value 
of a) a, there is residual degeneracy equal to the number of 
the magnetic flux quanta Nq. This degeneracy is resolved 
by a second pair of operators b, that commute with a, 
and depend on the guiding center coordinates of the electron, 
R a = r a — ^-j--K b i 2 B . Operators create the (unnormalized) 
single particle eigenstates of the lowest LL, 

C ■ M W = ^e- z * z/4 ^, (2) 

with z — x + iy being the complex coordinate of an electron 
in the plane (and z* denoting its complex-conjugate). The 
quantum number I is an eigenvalue of the angular momentum 
L z and the single particle states </>; are localized on concentric 
rings around the origin. 

To illustrate the effect of mass anisotropy, we take the prin- 
cipal axes of the mass tensor to be along the x and y directions 
(4> = 0), with different masses along the two directions (a ^ 
1). Via simple rescaling x — > x/y/a, y — > y\fa., and therefore 
introducing a oc \fom x + -^^y an d & oc y/omx ~ ^ 7r y we 
can immediatealy write down the single particle orbitals for 
this case: 

4>f{x >V )= [^=+iy^j e~^ 4+ay2 \ (3) 

Notice that the probability density \<pi | 2 is no longer localized 
on a circle, but rather an ellipse for a ^ 1, Therefore, on a 
single particle level, the effect of mass anisotropy is to stretch 
or squeeze the one-body orbitals along certain directions, pos- 
sibly rotating the principal axis (for <fi ^ 0). 

As we mentioned in Sec.Q] certain semiconductor materials 
are likely to have non-trivial metric defined by the anisotropy. 
Alternatively, the effective mass tensor can be experimentally 
tuned by tilting the magnetic field^. Tilting is known to 
produce complicated effects because it induces the coupling 
between electronic subbands and LL mixing, and a detailed 
analysis will be presented elsewhere. However, with a nor- 
mal confinement (perpendicular to the Hall surface) given by 
a harmonic well, the Hamiltonian with a strong perpendicular 
magnetic field B and a non-zero in-plane magnetic field Bn 
can be solved exactly, yielding two characteristic harmonic 
oscillator frequencies^: 

where uj z is the harmonic frequency of the normal confine- 
ment. The mixing of the cyclotron frequency and the confine- 
ment frequency is parametrized by tan 26 = 2ui c u^/ \u>1 — 
lo 2 z ), where ojii = eBu/m. Writing A, = i,2 = u)i/oj c , 
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the anisotropy parameter becomes a = AiA2/(Ai sin 2 9 + 
A2 cos 2 9) . Parametrizing the in-plane magnetic field by B y = 
(B\\ cos 0, sin 0), the effective metric associated with the 
tilt is given by 



cosh 29 + sinh 29 cos 2(f> 
sinh 2$ sin 20 



sinh 29 sin 2(f) 
cosh 29 — sinh 29 cos 2<\> 



(5) 



where cosh 29 = \(a+ Therefore, the effect of tilting on 
the LLL single-particle levels can be captured by the variation 
of the mass tensor. 

In order to study a finite, interacting system of N e elec- 
trons, it is convenient to choose a compact surface to rep- 
resent the 2DES. As we emphasized in Sec|I] the presence 
of mass anisotropy destroys rotational invariance, and one 
must use periodic boundary condition a 23,24 i.e. put the 2DES 
on the surface of a torus. The unit cell can generally be 
chosen as a parallelogram with sides a and b whose area 
S is quantized because of the magnetic translations algebra: 
S = I a x b| = 2ir£ 2 B N ( p. The single-particle states compati- 
ble with periodic boundary conditions are given in the Landau 
gauge by 

<f>f n (r) = ^— ^2 c i(x T +ka)v ~^ (Xi+ka+x)2 



H n | ^ + t a + X ) , (6) 



where j = 0, . . . , — 1, Xj = 2irj/b, normalization factor 
is Af = b^pK \fa2 n n\, and the sum over k extends over all 
integers. We have set £b = 1. The wavefunction for the nth 
LL involves a Hermite polynomial H n . For simplicity, we 
assumed the case of a rectangular torus and g = diag [a, 1/ct], 
but Eq.© can be generalized to an arbitrary shape/anisotropy 
using Jacobi theta functions. 

Many-body states, like in the isotropic case, can be classi- 
fied using a crystal quasimomentum KM defined in a Brillouin 
zone. With a suitable definition of the Brillouin zone, incom- 
pressible states always occur at K = 0, and are characterized 
by the gap in their excitation spectrum. The Hamiltonian for 
N e electrons is given by the sum of the kinetic term and the 
Coulomb interaction, 



H 



1 



(7) 



2DES is embedded in a three-dimensional dielectric host ma- 
terial, which is characterized by its own dielectic tensor e, 
defining the metric for the distance |r, — r_j| e between interact- 
ing electrons. This tensor in general can be different from the 
one that parametrizes the cyclotron orbits (a,<f>). However, 
physical properties are determined only by the relative differ- 
ence between the mass tensor and the dielectric tensor, and for 
simplicity we can set the latter to unity. In other words, we as- 
sume that Coulomb interaction is isotropic in space, hence its 

Fourier transform is V(q) = 1/q = 1/wg 2 + (A (to model 



finite-width effects, we use the softened form of V(q), fol- 
lowing the Fang-Howard prescription). Projected to a sin- 
gle nth LL, the interaction part of the Hamiltonian becomes 

h = y.\j : VhhhiAA^ ^ where 



v, 



^£V<q)£„(^) 



C l%,' ,! iX ' V »)<£ ■ ■ 6 , • ■ 

^J1-J4 U J1+J2,J3+J4' 



(8) 



where g 2 = g ab q a qb (e.g. in case of a diagonal mass tensor, 
g 2 = ag 2 + g 2 /a) and C n is the Laguerre polynomial. The 
primed <5-functions are to be taken (mod N^), and the sum 
over q extends over the reciprocal space (the prime on the 
sum indicates that the diverging q = term is cancelled by 
the positive background charge). 

Apart from the many-body translational symmetry, discrete 
symmetries can be used to further reduce the Hilbert space. 
Several types of Bravais lattices are possible, depending on 
the angle 9 between the sides of the torus, a and b, and the 
aspect ratio, |a|/|b|. Both of these can be tuned as free pa- 
rameters. In the presence of anisotropy, however, the highest 
symmetry is only given by the rectangular lattice, even when 
|a| = |b|. Tuning the angle 9 enables to perform the area- 
preserving deformations of the torus, which is useful in re- 
solving the collective modes of FQH states in finite systems, 
and probing quantities such as Hall viscosity^ 2 -. 

In this paper we only consider spin polarized electrons 
and neglect the so-called multicomponent degrees of freedom, 
which can be the usual spin or bilayer/valley degree of free- 
dom. This means that the filling factors we refer to as n + v 
correspond to kn + v in experiments, where integer k denotes 
the additional degeneracy that comes from several "flavors" of 
electrons. 



III. ANISOTROPY IN THE LOWEST LANDAU LEVEL : 
ROBUSTNESS AND THE INTRINSIC METRIC OF THE 
LAUGHLIN STATE 

In Fig.Q] we present the energy spectrum of the Coulomb 
interaction at v = 1/3 as a function of anisotropy (we assume 
4> = 0). The system is placed on the torus with a square unit 
cell, and energies are expressed in units of e 2 /eis- A very flat 
minimum around isotropy point and the existence of a robust 
gap suggest that the ground state of the Coulomb interaction 
at v = 1/3 is remarkably stable to variation in anisotropy. 
As we show below, in this range of a, the ground state is de- 
scribed by a generalized Laughlin wavefunction. Moreover, 
the set of lowest neutral excitations, forming a magneto-roton 
branch, are also stable and separated from the rest of the spec- 
trum. Within this manifold, some level crossings occur as a is 
changed, but this only corresponds to the redistribution of the 
levels within a roton branch. Beyond a w 0.5, the ground- 
state energy rises, indicating an instability and the eventual 
destruction of the Laughlin phase. 

In rotationally-invariant situations, the incompressible liq- 
uids at fillings ^=l/mofn = LL (m being an odd integer) 
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FIG. 1. (Color online). Energy spectrum in units of e 2 /ds as a 
function of anisotropy a for the square unit cell and n = LL 
Coulomb interaction at v = 1 /3. The system is iV e = 7 electrons 
and (j> = 0. Due to the square unit cell, the spectrum is symmetric 
under a — 1/a. 



are described by the Laughlin wavefunction 2 . In the geometry 
of an infinite plane, the Laughlin state is given by 



v — 1/m 



IK' 

i<j 



■ J2k z k z kl 



(9) 



Here z = x + iy stands for the usual complex representa- 
tion of the coordinates in the plane, and can also be expressed 
in terms of spinor coordinates providing a mapping to the 
spherical geometry^ 2 -. However, the Laughlin state can also 
be extended to the torus geometry^, where continuous rota- 
tion symmetry is broken down and survives at most in form 
of a discrete subgroup. In this case, ^> v ~ 1 l m is defined by 
its short-distance correlations which assume the form of the 
(odd) Jacobi -d\ theta function of — z^ 3 -. 

As a trial wavefunction, vf^ -1 ' 3 provides an excellent de- 
scription of the physical system at v = 1/3 in the limit of 
strong cyclotron energy with respect to the Coulomb repul- 
sion in the 2DES plane, haj c ^> e 2 /e£s, when excitations to 
higher LLs are prohibited. In this limit, the operators a, a act 
trivially and the only dynamical degrees of freedom are the 
(non-commuting) guiding centers. Up to the normalization, 
we can then view the wavefunction ((9) as follows 



*r 1/m (9)«n (bl(g)-b](g)) |0> 



(10) 



where b\(g) explicitly depend on the metric g 32 . For gen- 
eral g, b\ (g) is obtained by a Bogoliubov transformation from 
the bi, b\ in the rotationally invariant case. Equivalently, the 
wavefunction can be expressed by a unitary transformation 
^h{g) = exp(— i^ a ^A a/3 )^i(0), where £, a p is a real sym- 
metric tensor and A Q/3 = j ^2i{Rf, -Rf} is the generator of 
area-preserving diffeomorphisms 32 . The expression for the 
transformation matrix and the first-quantized expression for 
the wavefunction in Eq.dTOb is given in Ref.l28l 
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FIG. 2. (Color online). Overlap between the Coulomb ground state 
at v = 1/3 for fixed anisotropy ao = 2, (f>o — and the family of 
Laughlin states parametrized by varying a, <j>. The system is N e — 9 
electrons on a hexagonal torus. 



The freedom in choosing g implies that the usual 
rotationally-symmetric Laughlin wavefunction is a represen- 
tative of a class of wavefunctions. However, being a topolog- 
ical phase, the physics of the Laughlin state does not depend 
on any given metric or lengthscales. Various wavefunctions 
(ff) differ from one another microscopically in terms 
of the shapes of their elementary droplets. For a FQH state 
at filling v = p/q (p, q are not necessarily co-prime), an el- 
ementary droplet is a unit of fluid containing p particles in 
an area that encloses q flux quanta. The incompressible state 
is a condensate of such elementary droplets. For example, 
at v = 1/3 we have a single particle occupying each three 
consequtive orbitals and preventing more particles from pop- 
ulating this region. In the language of root partitions and the 
Jack polynomials 34 , the Laughlin v = 1/3 state is defined by 
a root pattern 100100100100100 . . ., and therefore its elemen- 
tary droplet is 100. Note that these simple patterns only serve 
as labels for correlated wavefunctions that cannot be thought 
of as a simple crystal of electrons pinned at each third orbital 
and repelling each other via electrostatic forces. 

For a model wavefunction, the "gauge" freedom in b, W im- 
plies that the shape of elementary droplets changes with vary- 
ing g, but the basic physical properties remain invariant unless 
the anisotropy magnitude a becomes too large or too small. If 
a is such that the maximum effective separation between elec- 
trons along some direction is of the order of Ib or smaller, the 
FQH liquid correlations are expected to break down and CDW 
states might be favored (we will present some examples of this 
in Sec. lIVb . On the other hand, a generic state involves a com- 
petition of two different metrics - the cyclotron metric, g m , 
and the Coulomb metric gc, therefore it is best approximated 
by the Laughlin state with an intrinsic metric g that is gener- 
ally different from both g m and gc- Intrinsic metric g is the 
one that minimizes the variational energy E g , 



E g = 



(* L (g)\H(g m ,g c )\* L (g)) 



(11) 



(*l(s)I*l( 5 )> 
To find the intrinsic metric in the microscopic calculation, 
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FIG. 3. Dependence of the intrinsic metric a on the mass metric a 
(Coulomb metric is set to identity). 



we use a slightly different criterion that g should maximize 
the overlap with the exact ground state of H(g m , gc)- In the 
rotationally-invariant case, the ground state of the Coulomb 
interaction at v = 1/3 is known to have a remarkably high 
overlap with the Laughlin wavefunction. The overlap is de- 
fined as a scalar product between two normalized vectors, and 
in this particular case it is typically greater than 97%. There- 
fore, we expect the intrisic metric chosen to maximize the 
overlap also to minimize the correlation energy ( fTTT >. To ob- 
tain the anisotropic Laughlin states, we perform exact diago- 
nalization of the "V\ Hamiltonian" on the torus. This Hamil- 
tonian gives vf^ -1 / 3 as a unique and densest zero-energy 
ground state 21 . Note that any translationally-invariant inter- 
action can be expanded in terms of the Laguerre polynomi- 
als, V(q) = Y^m Kn£ m (q 2 ), where the coefficients V m are 
the Haldane pseudopotentials 2 ^. Truncating this expansion 
at the first term, Vi£i(q 2 ), singles out the strongest ("hard- 
core") component, which defines the Laughlin state at filling 
v = 1/3. As the pseudopotential Hamiltonian is just a pro- 
jection operator in the relative angular momentum space, the 
metric in ViC\(o^) is the same as that originating from the 
cyclotron orbits. 

In Fig. |2] we pick the ground state of the Coulomb in- 
teraction with fixed mass anisotropy chq = 2,<po = (the 
metric of the dielectric tensor is implicitly assumed to be 
a = l,(f> = 0), and we evaluate the overlap with a family 
of Laughlin states generated by varying a, <j>. The overlap 
|(^|^ O=2 ' 0O=O )| is plotted as a function of a and 0. We 
observe that the principal axis of the Laughlin state is aligned 
with that of the Coulomb state (maximum overlap occurs for 
<j> = <pQ = 0). Interestingly, the maximum overlap does not 
occur for a = ao, but for some value of the anisotropy that is 
a "compromise" between the dielectric a = 1 and a cyclotron 
one a — 2. The value of the anisotropy that defines the in- 
trinsic metric depends linearly on the band mass anisotropy 
(Fig- El). This result illustrates the ability of the Laughlin state 
to optimize the shape of its fundamental droplets and maxi- 
mize the overlap with a given anisotropic ground state of a 
finite system. 

An alternative way to obtain the intrinsic metric is to an- 
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FIG. 4. (Color online). Energy spectrum of N e = 9 electrons at 
v — 1/3 with the effective mass anisotropy a m = 2 along the x- 
axis. When plotted as a function of \/kT+ky (green crosses), two 
branches of the magneto-roton mode are present (blue dotted lines 
are guide to the eye). If the spectrum is plotted as a function of 
\J g ab k a kb, the two branches collapse onto the same curve. 



alyze the shape of the lowest excitation - the magneto-roton 
mode, which was successfully described by the single-mode 
approximation 3 ^. In a rotationally-invariant case, this mode 
has a minimum at k ~ £~ . In the presence of anisotropy, 
the minimum occurs at different | k | in the different directions 
(Fig. 2|. This leads to an alternative definition of the intrin- 
sic metric based on the shape of the roton minimum in the 
2D momentum plane. We numerically establish that this defi- 
nition agrees well with our previous definition of the intrinsic 
metric. In Fig.|4]we plot the energy spectrum of an anisotropic 
Coulomb interaction at v = 1/3 as a function of the rescaled 
momentum \J g ab k a k^ where g is the guiding center met- 
ric that maximizes the overlap with the family of Laughlin 
wavefunctions (Fig. |3)- With the usual definition of the mo- 
mentum |fc|, several roton minima appear. Different magneto- 
roton branches collapse onto the same curve if we plot them 
as a function of \J g ab k a kb. This is reasonable, because the 
magneto-roton mode is well approximated by single-mode 
approximation up to the roton minimum 36 , which is defined 
entirely in terms of the properties of the ground state. The 
anisotropy of the ground state structure factor (determined by 
the shape of elementary droplets) dictates the position of the 
roton minimum. 

The analysis of this section in principle applies to other 
LLL states at fillings v = pj (2p + l),p = 2, 3, . . ., though it 
is more involved because of the "multicomponent nature" of 
these states and typically a smaller excitation gap. 



IV. HIGHER LANDAU LEVELS : QUANTUM PHASE 
TRANSITIONS DRIVEN BY ANISOTROPY 

We have found that v = 1/3 in the LLL is particularly ro- 
bust with respect to anisotropy, and this is the case also with 
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FIG. 5. (Color online). Energy spectrum of v = 1/3 (left) and 
v — 1/2 filled n — 2 LL (right): mass anisotropy establishes and 
reinforces the stripe order. 



other prominent FQH states. In higher LLs, due to a number 
of nodes in the single-particle wavefunction, the region of the 
phase diagram where incompressible states occur becomes in- 
creasingly narrower, and compressible phases such as stripes 
and bubbles take over. In this Section we discuss the effects of 
anisotropy on FQH states in higher LLs, focusing on fillings 
v = 1/3 and 1/2. Because of closer energy scales, we find 
that moderate changes in the anisotropy induce phase transi- 
tions between compressible and incompressible phases. 



A. n > 2 Landau levels: stripes and bubbles 

In n = 2 LL and higher, isotropic FQH states are energet- 
ically less favorable than stripe and bubble phases at filling 
v = n + 1/2 and v = n + 1/3, respectively. In Fig. [5] we 
show the energy spectrum (in units of e 2 /els) as a function 
of anisotropy a (we set the angle <f> to zero). Energies are 
plotted relative to the ground state at each a, and we chose 
the relatively modest system sizes (N e = 8 and 10 electrons) 
to facilitate comparison with the existing isotropic data in the 
literatur e) 16 - 17 . The aspect ratio is set to the optimal values for 
the appearance of stripes or bubbles (see Refs.[l6land[l7l). 

As we see on the right panel of Fig. [5] at v = 1/2 the 
presence of mass anisotropy reinforces the stripe when a is in- 
creased. This leads to a more pronounced quasi-degeneracy of 
the ground-state multiplet, and an increase of the gap between 
this multiplet and the excited states. For yet larger values of 
a, it appears that some of these excited states may become 
the ground state, however this occurs for very large a when 
this finite system effectively becomes one-dimensional under 
anisotropy deformations. 

In case of v = 2 + 1/3 state, it has been argued that the 
isotropy point is described by a two-dimensional CDW order 
known as the bubble phase 17 . A bubble differs from a stripe 
in having a larger degeneracy and a two-dimensional mesh 
of (quasi)degenerate ground-state wavevectors (as opposed to 
the one-dimensional array in case of a stripe). The spread 
of the quasidegenerate levels was also found to be somewhat 
larger than in case of stripes. All of these features are obvious 
in Fig |5] (left) for a = 1. The bubble phase remains stable to 
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FIG. 6. (Color online). Spectrum of N e — 8 electrons at v — 1 + 1/3 
with thickness w = 21b ■ Inset: same spectrum plotted relative to the 
ground state at each a. Unit cell has a rectangular shape with aspect 
ratio 3/4. 

some extent when a is reduced; for very small a it is even- 
tually destroyed and replaced by a simple CDW. On the other 
hand, when a is increased, a smaller subset of momenta be- 
comes very closely degenerate with some of the excited lev- 
els. This second-order (or weakly first order) transition re- 
sults in a stripe phase. As for the v = 1/2 case, this stripe 
becomes enhanced as a is further increased. Therefore, in 
n > 2 LLs mass anisotropy generally produces stripes, even 
when isotropic ground states have a tendency to forming a 
bubble phase. 

B. n = 1 Landau level: incompressible to compressible 
transitions driven by anisotropy 

In n = 1 LL, v = 1/3 state is significantly weaker than 
its n = LL counterpart, having an experimental gap an or- 
der of magnitude smaller and roughly the same as the gap of 
v = 1/2 state. This has been anticipated in early numerical 
calculations that found the ground state of the Coulomb in- 
teraction projected to n = 1 LL to be at the transition point 
between compressible and incompressible phases^. 

Although idealized numerical calculations with pure (pro- 
jected) Coulomb interaction work exceedingly well in n = 
LL, more realistic models are required to describe phases in 
?i = 1 LL. In particular, the inclusion of finite width effects^ 
and varying a few strongest Haldane pseudopotentials in nec- 
essary to determine the phase diagram. We find that varying 
the Vi pseudopotential leads to the following outcomes: (i) 
generically, for 5 Vi < 0, the system is pushed deeper into a 
compressible phase; (ii) for SVi > 0, finite-size calculations 
on systems up to N e = 9 electrons permit the existence of two 
regimes: for < SVf < 8V\ < 5V±, the ground state is in 
the Laughlin universality class, but the lowest excitation is not 
the magneto-roton; for 8V\ > SVi, the ground state and the 
excitation spectrum is the same as in n = LL. For smaller 
systems, SV^ is estimated to be around 0.1e 2 /e£s, while SV" 
is around 0.04e 2 / c£b- Larger systems suggest that these two 
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FIG. 7. (Color online). Guiding-center structure factor So(q) for 
v = 1/3 state in n = 1 LL with thickness w = 21b and anisotropy 
a = 0.4 (a). For comparison, we also show So(q) for the state with 
a = 1.3 which is in the Laughlin universality class (b). Two peaks in 
the response function (a) represent the onset of compressibility and 
CDW ordering. 



points might merge in the thermodynamic limit, when only a 
small modification of the interaction might be needed for the 
Laughlin physics to appear at^ = l/3inn=l LL. Alterna- 
tively, we can consider the Fang-Howard ansatz that mimicks 
the finite-width effects. In this case, the width of is or smaller 
is sufficient to drive a phase transition between the compress- 
ible state and the Laughlin-like state, in agreement with results 
on the sphere and using an alternative finite-width ansatz^. 

In summary, the ground state at v = 1 + 1/3 very likely 
belongs to the Laughlin universality class. We note that the 
collective mode in this case displays significantly more wig- 
gles than in the LLL (some wiggles exist in case of n = LL 
Coulomb state, but they are less pronounced). For large mo- 
menta, the magneto-roton mode also appears to merge with 
the continuum of quasiparticle-quasihole excitations. This is 
likely a finite-size artefact, although we cannot rule out that it 
represents an intrinsic feature of f = 1 + 1/3 state, in which 
case it might have an observable signature in optical experi- 
ments that distinguishes it from v = 1/3 state. 

Because of the fragility of v = 1 + 1/3 state, we expect that 
mass anisotropy might have more dramatic consequences than 
in the LLL. In Fig. [6] we plot the energy spectrum as a function 
of anisotropy. One notices that the isotropy point (a = 1) does 
not bear any special importance - indeed, the system appears 
more stable in the vicinity of it where it can lower its ground 
state energy or increase the neutral gap. On either side of the 
isotropy point, however, the system remains in the Laughlin 
universality class; e.g. at a = 0.8 and a = 1.3 the maximum 
overlap with the Laughlin state is 75% and 80%, respectively 
(these overlaps, although modest compared to the standards of 
n = LL, can be adiabatically further increased by tuning the 
V\ pseudopotential). Note that the quoted maximum overlaps 
are achieved by the Laughlin state with a' somewhat different 
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FIG. 8. (Color online). Spectrum of N e — 14 electrons at v = 
1 + 1/2 with thickness w — 21b, as a function of anisotropy a. 
Energies are plotted relative to the ground state at each a, and the 
unit cell has a rectangular shape with aspect ratio 3/4. 

from a of the Coulomb state, analogous to Fig|2j 

The new aspect of Fig|6]is the transition to a compressible 
state with CDW ordering for a < 0.4. In that region of pa- 
rameter space, the system is very sensitive to changes in the 
boundary condition - the sharp degeneracies seen in rectangu- 
lar geometry in Fig|6]are not obvious in case of higher sym- 
metry, square or hexagonal, unit cell. As an additional diag- 
nostic tool for the compressible states, it is useful to consider 
a guiding-center structure factor, 

Sb(q) = -^^(e jq - R 'e- tq - R ^)-(e iq - Rt )(e- iq ' Rj ), (12) 

^ 1,3 

where the expression for the Fourier components of the 
guiding-center density, p(q) = e* qIl \ has been used. 
Note that So(q) is normalized per flux quantum rather than 
(conventional) per particle^. In FigfTJa) we show the plot of 
So(q) evaluated for the state with a = 0.4 in Fig|6] Two 
sharp peaks in the response, similar to those previously iden- 
tified in n > 2 LL states^, are the hallmark of CDW order. 
They are to be contrasted with the smooth response in case 
of an anisotropic state in the Laughlin universality class for 
a = 1.3,Fig|3b). 

As a second example in n = 1 LL, we consider half fill- 
ing where the Moore-Read Pfaffian state 5 is believed to be 
realized in some regions of the phase diagram. This state 
has a non-Abelian nature, which is reflected in the non- 
trivial ground state degeneracy^ when subjected to periodic 
boundary conditions. For v = 1/2, the eigenstates of any 
translationally-invariant interaction possess a twofold center- 
of-mass degeneracy^. On top of this, Moore-Read state 
has an additional threefold degeneracy. Conventionally, the 
many-body Brillouin zone is defined for p = l,q = 2 and 
has a size N 2 (N being the GCD of N e and N<j,), which 
forces the degenerate groundstates to belong to a Brillouin 
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zone corner K = (N/2, N/2) and centers of the sides, 
K = (0, N/2); {N/2, 0). It is also possible to define a "quar- 
tered" Brillouin zone such that the three degenerate states are 
all mapped to zero momentum 38 . The three K sectors are 
equivalent for a hexagonal unit cell, however in an anisotropic 
system the degeneracy is always lifted. 

In FigJH]we plot the spectrum of the Coulomb interaction 
as a function of anisotropy (states belonging to K sectors 
where the Moore-Read state is realized, are indicated). As 
earlier, we assume finite width of w = 2£b in order to in- 
state the Pfaffian correlations^. Note that our calculation 
only uses two-body (Coulomb) interaction, therefore in each 
finite system the Moore-Read state will mix with its particle- 
hole conjugate pair, the anti-Pfaffian 3 ^. The mixing between 
the two states can be controlled by including higher LLs^S. 
For 0.5 < a < 1.3, we find a three-fold quasi-degenerate 
multiplet, suggesting the presence of Moore-Read state at the 
isotropy point and in the neighborhood of it. In finite systems, 
there is some splitting of the degeneracy that might be reduced 
upon tuning the Vi, V3 pseudopotentials. Also, upon tuning 
the anisotropy around a = 1, there are crossings within the 
multiplet of degenerate ground states without apparent clos- 
ing of the gap. The region of the Moore-Read state is defined 
by sharp transitions towards crystal phases. These transitions 
are likely second order because they do not appear to involve 
any level crossing, but rather lifting of the degeneracy within 
a ground-state multiplet. 

V. CONCLUSION 

We have presented a method to study the effects of 
anisotropy on FQH phases in finite-size systems. We found 
that the prominent FQH states (as in the lowest Landau level) 
are robust to variations of anisotropy, due to the adjustment 



of the intrinsic metric describing the shape of their elemen- 
tary droplets. As we demonstrated using the example of the 
Laughlin v = 1/3 state, this metric is usually a compromise 
between the metric dictated by the cyclotron motion and the 
metric originating in the dielectic environment of the 2DES. 
In this sense, it is unlike the non-interacting Landau level 
problem, or the problem of weak localization— where the 
anisotropy can be completely "gauged away" (i.e. removed) 
by length rescaling. Instead, it is more akin to the problem 
of shallow donors in many-valley semiconductors^ 3 -. Indeed, 
such compromise picture leads to a quantitatively accurate de- 
scription of the variation of the critical density of the metal- 
insulator transition (an intrinsically many-body phenomenon) 
in three-dimensional doped many -bailey semiconductors^, so 
one may wish to ascertain to what extent this can lead to 
quantitative predictions in the FQHE case. In higher LLs, 
anisotropy induces quantum phase transitions, likely of sec- 
ond order, to compressible phases with broken symmetry. 

Anisotropy is an important aspect of FQHE as it represents 
a mechanism that probes the intrinsic metric of incompress- 
ible fluids in the geometrical picture of the FQHE. In addi- 
tion, because our calculations show the possibility of phase 
transitions in the n > Landau levels as a function of mass 
anisotropy, it motivates experimental studies on systems with 
both moderate mass anisotropy (e.g. AlAs and Si, a ~ 3 — 5), 
as well as systems with large mass anisotropy (e.g. Ge, 
a ~ 20), where behavior may be different in the upper Lan- 
dau levels from the anisotropic GaAs. In these systems, as 
in GaAs, anisotropy could be furher tuned using tilted fields, 
thereby adding to the richness of the FQHE phenomena. 
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